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Abstract 


In this paper we have proposed an almost unbiased estimator using known value of some 
population parameter(s). Various existing estimators are shown particular members of the 
proposed estimator. Under simple random sampling without replacement (SRSWOR) scheme the 
expressions for bias and mean square error (MSE) are derived. The study is extended to the two 
phase sampling. Empirical study is carried out to demonstrate the superiority of the proposed 
estimator. 


Key words: Auxiliary information, bias, mean square error, unbiased estimator, 
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1. Introduction 


Consider a finite population U = U41, Us, ...Un of N units. Let y and x stand for the 
variable under study and auxiliary variable respectively. Let (yi, xi), i21, 2,.., n 
denote the values of the units included in a sample s, of size n drawn by simple 
random sampling without replacement (SRSWOR). The auxiliary information has 
been used in improving the precision of the estimate of a parameter (See 
Cochran (1977), Sukhatme et. al. (1984) and the references cited there in). Out 
of many methods, ratio and product methods of estimation are good illustrations 
in this context. 


In order to have a survey estimate of the population mean Y of the study 
character y, assuming the knowledge of the population mean X of the auxiliary 
character x, the well-known ratio estimator is 


(1.1) 


Bahl and Tuteja (1991) suggested an exponential ratio type estimator as — 


[a = V CX X-x 
re y Xa 





Pak.j.stat.oper.res. Vol.IV No.2 2008 pp63-76 


Rajesh Singh, Mukesh Kumar, Florentin Smarandache 


several authors have used prior value of certain population parameter(s) to find 
more precise estimates. Sisodiya and Dwivedi (1981), Sen (1978) and 
Upadhyaya and Singh (1984) used the known coefficient of variation (CV) of the 
auxiliary character for estimating population mean of a study character in ratio 
method of estimation. The use of prior value of coefficient of kurtosis in 
estimating the population variance of study character y was first made by Singh 
et. al. (1973). Later used by Singh and Kakaran (1993) in the estimation of 
population mean of study character. Singh and Tailor (2003) proposed a modified 
ratio estimator by using the known value of correlation coefficient. Kadilar and 
Cingi (2006), Khosnevisan et. al. (2007), Singh et. al. (2007) Singh and Kumar 
(2009) and Singh et. al. (2009) have suggested modified ratio estimators by 
using different pairs of known value of population parameter(s). 


In this paper under SRSWOR, we have proposed almost unbiased estimator for 
estimating Y. 
2. Almost unbiased ratio type estimator 


Suppose 
| _ | — (aX+b —  — ((aX+b)-(ax+b 
to — Us = »( ) Use = exp lere) 





x+b 


ouch that ty t,,,t,..€ Ww, where w, denotes the set of all possible ratio type 
estimators for estimating the population mean Y. By definition the set w, is a 
linear variety, if 

twr = OY t Gt. t Qt. € W (2.1) 
for Z,a 21, ER (2.2) 
where w; (i20, 1, 2) denotes the statistical constants and R denotes the set of 
real numbers. 


To obtain the bias and MSE of t, we write 
y=Y(l+e,), x=X(l+e,), 











such that 
E (eo) =E (e1)-0. 
E(e9) = fiC;. E(e?) 2 fqC?, E(ege,) = fipC, C, 

1 1 1 oq me 1 4 me 
here f, = (——- —), B e -Y), S= (-XJ, 
where fj =-=) (x 240^ | woe | 

S S C S 
ap a C= K=p| — |, p=, 
yy ^ x T3 " Ss.) 
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Expressing t in terms of e’s, we have 


ty =¥(1 + e0) [wo + ,(1 + 0e,)"1 + wexp f- 5: (1+ 86,)71]] (2.3) 
aX 
where 8 = 
aX+b 


Expanding the right hand side of (2.3) and retaining terms up to second order of 
e's, we have 


t, =Y E +e, — «8e, + 67 (o t =) ei — Bwege] (2.4) 
where w= (w, +=). 

Taking expectations of both side of (2.4) and then subtracting Y from both side, 
we get the bias of the estimator tẹ, up to the first order of approximation as 








B(t,) — fY [a?c2 (o. es) = BwpC, C,| (2.6) 
B(t,,) = f, Y[e?c? — pc, c,] (2.7) 
30“ 5 apC, C. 
B(t ER f,Y = =] (2.8) 
From (2.4), we have 
(t, — Y) = Y[e, — Owe, | (2.9) 


Squaring both sides of (2.9) and then taking expectations, we get MSE of the 
estimator t, up to the first order of approximation, as 


MSE(t,,) = f,Y [CZ + €? C2 — 280pC, C] (2.10) 
This is minimum when 

o — k(- p) (2.11) 
Putting this value of (= k) in (2.10), we get the minimum MSE of t,, as 

min. MSE(t,,) = f, Y?C2(1— p?) (2.12) 


which is same as that of traditional linear regression estimator 
from (2.5) and (2.11), we have 


(4 = | (2.13) 


From (2.2) and (2.13), we have only two equations in three unknowns. It is not 
possible to find the unique values for w; i = 0,1,2. In order to get unique values 


for «»,.,, we shall impose the linear restriction 
cogB(y) + w,B(t,,.) + o,B(t.,) = 0 (2.14) 
a (2. E^ (2.11) and (2.14) can be written as in the matrix form as 
1 1 
b a i5 Is ji (2.15) 
0 B(t,.) Bit...) 
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Using (2.15) ,we get unique values of œw, (i=0,1,2) as 

Ag 

Ay 
- 
where 

A,— Bt...) - B(t,.) 
_ poe wai ous fu 

Aro” B(t,.. {1 — k} + 2 B(t,..) fk | | (2.1 7) 


A= k.B(t,,..) 
A,,— —k B(t,,) 


Use of these w;, (i=0,1,2) remove the bias up to terms of order (n'') at (2.1). 


3. Product -type estimators 


e _ ( (a¥+b)—(aX+b) 
SUB =y tose = aed 





such that to, the, tose €Q, where Q denotes the set of all possible product -type 
estimators for estimating the population mean Y. By definition, the set Q is linear 
variety if 


ta = doy T qitps t Gatyce EQ (3.1) 
for 32.,4;=1, QER (3.2) 


where q,(i=0,1,2) denotes the statistical constants. 


Expressing t, in terms of e’s, we have 


ta = $(1- ej) [oo e (14 864) + o; exp (Ps (1 + Be,)~}| (3.3) 
where 6 = =, 
aX+b 


Expanding the right hand side of (3.3) and retaining terms up to second power of 
e's, we have 


t, = Y|1-- ej +8qe, — e? + q&e,e, | (3.4) 
where q = q; + 2 (3.5) 
Taking expectations of both sides of (3.4) and then subtracting Y from both 


sides, we get the bias of the estimator t,, up to the first order of approximation 
as 


B(t,) = f,Y [-2 8C? + q8pC, C " (3.6) 
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Bias expression for the estimators t,, and  t,,,, is given by 
B(t,.) =f, Y [@pC,c, | (3.7) 

ES l 42 Bp yy 
B(t s) = f,Y |-:6 c? H] (3.8) 
From (3.4), we have 
(t4 — Y) = Yle + 8gei] (3.9) 


Squaring both the sides of (3.9) and then taking expectations, we get MSE of the 
estimator t,, up to the first order of approximation, as 


MSE(t, ) = £Y? [c2 + 62q2c2 + 26qpC, C, ] (3.10) 
which is minimum for 

q= -k- -pd (8.11) 
Putting this value of q(=-k) in (3.10),we get the minimum MSE of t, as 

min. MSE(t,) = f, Y^C2(1 — p°) (3.12) 
which is same as that of traditional linear regression estimator. 


From (3.5) and (3.11), we have 
dE —— (3.13) 


From (3.2) and (3.13), we have only two equations in three unknowns. It is not 
possible to find the unique values for qs, i20,1,2. In order to get unique values of 
qi s, we shall impose the linear restriction 


qB (F) T qB(t,.) Li q; B(t,..) =0 (3.1 4) 


2 m (3.2),(3.13) and (3.14) can be written in the matrix form as 


1 
0 ui ^ df :H- (3.15) 
0 ss.) B(tpse) 


Solving (3.15), we get the unique values of qs (i20,1,2) as- 


Apo 
do = FR 
qı = = (3.16) 
_ Ap2 
q2 "" 
where 
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= B(t,..) - B(t,.) 
ASQ "o +k}+ B(t)[-k- 4} | 
Ani= —k BIER) 

Ap2= K B (tps) 


(3.17) 





Use of these q's (i=0,1,2) remove the bias up to terms of order o(n'') at (3.1). 


In Appendix A we have listed some of the important known estimators of the 
population mean, which can be obtained by suitable choice of constants 
0,,1—0,1,2, q,,i=0,1,2 and a and b. 


4. Proposed estimators in two phase sampling 


When X is unknown, it is sometimes estimated from a preliminary large sample of 


size n' on which only the characteristic x is measured (for details see Singh et. al. 
(2007)). Then a second phase sample of size n (n < n’) is drawn on which both y 


and x characteristics are measured. Let x — LX -,X, denote the sample mean of 
x based on first phase sample of size n', y = T -,y, and x= Ya be the 
sample means of y and x respectively based on second phase of size n. 

In two phase sampling the estimator t,,. will take the following form 

twa = Opa yt wid trasd T O»^q Used € Wa (4.1 ) 
for EZ. ua = 1, WER (4.2) 


(aX +b) (ax+b 
where tps = y(=*) and trse = yex exp [Ec 


To obtain the bias and MSE of t,,. we write 
y= Y(l+e,), x=X(l+e,), x 2 X(14 ej) 
such that 

E (eo) = E (e1) = E(e1) =0. 


E(e$) =f Cy, E(ej) -fjC$,  E(e'f)=fC}  E(eoei)- fipC, C, 
E(ege' ) = fopC Cx E(eje, ) — DC. 

1 1 1 1 
where fi A 1o CONI 


Following the procedure mentioned in section 2 and 3, we get bias and MSE of 
twa as 


B (toa) = y [e^ Cf; (wg + 8) — Op Cy C.F (044 + 225) (4.3) 
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MSE(t..4) = Y [ACZ + f; C293 — 2f eC, C, wa] (4.4) 

1 1 
where f; = -——=f, — fy. 

n n 
MSE (t,,4 ) is minimum, when 

wd 

Wig +- — 04 — k. (4.5) 
Putting this value of w, in (4.4), we get the minimum MSE of t,,; as 
min. MSE(t a ) = ¥°C2[f, —f,p7] (4.6) 
This is same as that of traditional two phase linear regression estimator. 
The bias expression for the estimators t,., and t,.4. is respectively given by 
B(t,4) = Y[6?C2f£,0,4 — 99 C, C. f, 0,4] (4.7) 
B(treae) = Y [e C2f, 792: — 8C, C, f, = (4.8) 


From (4.2) and (4.5), we have only two equations in three unknowns. It is not 
possible to find the unique values for wig’s i=0,1,2. 


In order to get unique values of w;,, we shall impose linear restriction 








Woa B(y) T @,4B(t,,4) T C) B(trsed) =0 (4.9) 
vanes (4.2), (4.5) and (4.9) can be written in matrix form, as 
1 
0 t A | m k! (4.10) 
0 a(t...) B(t,.e4) 
Solving (4.10), we get the unique values of w,.s, (i = 0,1,2) as 
Apa 
Ogg = — 
od Ara 
_ Asa | 
Oa = 5 (4.11) 
O24 = = 
ra 
where 
åa- B (trasd) 7 > B (trea) 
| | 1 ! 1 
Aroa” B(t sea) {1 7 k} + 2 B(t,.a) fk- =} 3 (4.12) 


Aua= k. B(t sad) 

Ava —k .B(t a) 
Use of these wie's (i=0,1,2) will remove the bias up to terms of order O(n") at 
(4.1). 
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The estimator t, written in (3.1), in two phase sampling, will take following form 
tog = Qoa + didtpsd + Azdtpsed € Qa (4.13) 
For Eod = 1, aa ER. (4.14) 
where q;,(i=0,1,2) denotes the statistical constants . 
The estimators t,.3 and t,.ge are 

_ f ax+b 
tosa = y) and 


t {== =| 
psed — Yexp (ax+b)+(ax'+ 


Following the procedure of section 4, we get the unique values of digs (i=0,1,2) 











as 

| Spod 

od " 

dia = yy (4.15) 
_ Ap2a 

2d "e 

where 


AST Btm) i : B(t..) 
Anoa= B(traea ){1 +k} + E B(t,.a) [7k B 3 | (4.16) 
Apid^ -k B(treea ) 
Anza =k.B (trea ) 


where 
B(t,.4) E Y[0a;.f;0C, C. | (4.17) 
B(tpsed ) = Y[8S5pc,c, — 7 qa fC?| (4.18) 


The minimum MSE of taa is given by 
MSE(t,4) = Y°C>[f, — f,p°) . 


5. Empirical study 


For empirical study we use the data sets earlier used by Kadilar and Cingi (2006) 
(population 1) and Khosnevisan et. al. (2007) (population 2) to verify the 
theoretical results. 
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Data statistics 
TTS 77 6 [sw 
Population? | 20 | 8 | 1965 | 188 | 0.956 | 0004 | 092 | 306. 





Table 5.1: Values of wis and qi's 


|o ws Population 1 Population 2 






8.590718 7.892148 
q ( 21.417) ( 2.919085) 
"m 11.86615 5.234461 
Ch (16.14158) ( 3.576773) 
"m -19.4569 -12.1266 
TE ( -36.5586) ( -5.49586) 
(PRE) 


The percent relative efficiencies (PRE) of various estimators of Y are computed 
and presented in Table 5.2 below. 


Table 5.2: PRE of different estimators of Y 


| Estimator | PRE (Pop I) | Estimator | PRE (Pop II) | 
— td | 10 | q | 100 | 
| t (| 212810 | qe | 4655001 — 
| t  — (| 143928 | œ | 285920 — 
| to | 143988 | œ | 338487 | 
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In order to see the performance of the suggested estimators in two phase 
sampling we use the data set of Murthy (1967) (Population Ill) and Steel and 
Torrie (1960) (Population IV). 


Population | C [| & |p | N- 


[Population 1 | 0.3542 | 0.9484 | 0.9150 | 80 | 30 | 10 | 
Population2 0.4803 | 0.7493 | -0.4996 


Table 5.3: The values of w;4's and q;,.'s 





Population Ill Population IV 
gg -0.241523 3.011435 
(doa) (1.808833) (1.089979) 
(4g -0.558071 1.370950 
(di) (0.125381) 0.730464) 


1.799595 -3.382385 
(-0.934214) (-0.820443 ) 





The percent relative efficiencies of various estimators of Y in two phase sampling 
are computed and presented in Table 5.4 below. 


Table 5.4 : PRE of different estimators of Y in two phase sampling 


| Estimator | PRE (Population I) | PRE (Population Il) 
36.642 24.562 
4.849 59.770 
200.420 48.365 

t 


23.628 115.142 
276.156 63.452 


123.762 
123.762 





6. Conclusion 


From theoretical discussion and empirical study we conclude that the proposed 
estimators under optimum conditions perform better than other estimators 
considered in the article. The relative efficiencies of various estimators are listed 
in Table 5.2 and 5.4. 
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Appendix A 


Table A.1: Some members of the proposed family of estimators - 


Ratio Estimator 


(corresponding to 
€); i=0,1 52) 


The mean per unit estimator 


A 
t,-—v- 
1757y- 


The usual ratio estimator 


x+C, 


Sisodia and Dwivedi (1981) 
estimator 


t „X +B) 


Singh et. al. (2004) estimator 
V XB2(x)+Cx 
xp» (x) + C, 
Upadhyaya and Singh (1999) 
estimator 
_— XC, * P260 
XC, +B (x) 


Upadhyaya and Singh (1999) 
estimator 


| _|X+p 
ts =F Zio 


Singh and Tailor (2003) 
estimator 


(7; =yex X-x 
pet d Erat 


Bahl and Tuteja (1991) 


estimator 
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Product Estimator 
(corresponding to q;.i=0,1,2) 


do — y 


The mean per unit estimator 
x 

qi~ yg 

The usual product estimator 

-— e 

C X«C, 


Pandey and Dubey (1988) 
estimator 


q2 


_X+B5(x) 
q3 => 
X +B (x) 

Singh et. al. (2004) estimator 
v xp» (X) "P Ce 

Xp»5(x) * C, 
Upadhyaya and Singh (1999) 
estimator 
qs = _ xC, + B(x) 

XC, s B(x) 


Upadhyaya and Singh (1999) 
estimator 


ET 


Singh and Tailor (2003) 
estimator 


— yex x-X 
q7 y P xax 


q4 


Bahl and Tuteja (1991) 


estimator 
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Ratio Estimator 


(corresponding to 
i=0,1,2) 


l 


pcm 
E "FEX x25, 0) 


Singh et. al. (2007) estimator 


to = vex _ X-x 
ind ERR. 


oingh et.al. (2007) estimator 


tin = Vexp| = 
107 Y dex. 


oingh et. al.(2007) estimator 


exp OK 
a1 BOX 
Singh et. al. (2007) estimator 


1 exp 9 
j CXV BA 


Singh et. al. (2007) estimator 


_ C (X-X) 
l =YeEXP ——= n 
C (X+x)+2p 


Singh et. al. (2007) estimator 


Singh et. al. (2007) estimator 


B. C0 (X —3) | 
B, C9 942p 
Singh et. al. (2007) estimator 


lis =Y E 
XX--x)2f, (x) 


Singh et. al. (2007) estimator 


[,;— yex 


Product Estimator 


(corresponding to q,,i=0,1,2) 


: edi io 
X+X+26, (x) 


Singh et. al. (2007) estimator 


ae Zr FFC. | 

oingh et.al. (2007) estimator 
exp AA 

dio = Yexp ene 2p 


oingh et. al.(2007) estimator 


BOK-X) 
BOX) 
oingh et. al. (2007) estimator 


l| Ga% 
Mad pes C (X+ +2300 


singh et. al. (2007) estimator 


-— C (x — X) 
da» P CX 342p 


q, j^Yex 


oingh et. al. (2007) estimator 


px - X) | 


Cx 
= yep EC 


oingh et. al. (2007) estimator 


8,2) @—X) | 
B.C) (X 922p 
oingh et. al. (2007) estimator 


- pa-X) 
dis YON TX 428,00) 


Singh et. al. (2007) estimator 


q,,—yex 
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In addition to above estimators a large number of estimators can also be 
generated from the proposed estimators just by putting different values of 
constants wi = 0,1,2, q;,i=0,1,2 and a and b. 
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